Magnetically induced electronic ferroelectricity in half-doped manganites 
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Using a joint approach of density functional theory and model calculations, we demonstrate 
that a prototypical charge ordered half-doped manganite, La^Cax^MnOa is multiferroic. The 
combination of a peculiar charge-orbital ordering and a tendency to form spin dimers breaks inversion 
symmetry, leads to a ferroelectric ground-state with a polarization up to several [iC/cm 2 . The 
presence of improper ferroelectricity does not depend on hotly debated structural details of this 
material: in the Zener-polaron structure we find a similar dramatic ferroelectric response with a 
large polarization of purely magnetic origin. 



PACS numbers: 71.15.Mb 75.47.Lx 77.80.-e 

Materials with simultaneous magnetic and ferroelec- 
tric ordering -multiferroics- are attracting enormous sci- 
entific interest [HE]- They offer the potential to control 
the magnetic orderparameter by the ferroelectric one and 
vice versa -a very desirable property from a technological 
point of view [3] . Such control requires large multiferroic 
couplings: a substantial ferroelectric polarization needs 
to be induced by the magnetic ordering. Even if in quite 
a few materials ferroelectricity and magnetism coexist, 
multiferroic couplings are tiny [H \5\ O E] . When design- 
ing materials with large multiferroic couplings, one has to 
exclude from the start the largest class of multiferroics, 
the ones in which multiferroicity is driven by spiral mag- 
netic ordering. In these materials multiferroicity relies on 
relativistic spin-orbit coupling as a driving force, which 
is intrinsically weak [T]. Charge ordered magnetic com- 
pounds are a far more promising class of materials with 
potentially large multiferroic couplings [8]. Coexistence 
of charge ordering and magnetism is found in a substan- 
tial number of transition metal oxides |9J , but to become 
strongly multiferroic a material needs to meet three ad- 
ditional requirements: (i) the symmetry is such that the 
magnetic ordering can push the charge ordering pattern 
from site-centered to bond-centered or vice versa [10] . {it) 
the material is insulating, as it has to support a ferroelec- 
tric polarization and (m) the material is electronically 
soft, so that inside it charge can easily be displaced. Half 
doped manganites of the type Lai/2Cai/2Mn03 famously 
meet the last two requirements Here we show they 
also fulfill the first one, allowing strong multiferroic cou- 
pling to emerge. 

In a combined approach of ab-initio density functional 
and model Hamiltonian calculations we show that the 
strong electron-electron interactions together with the 
Jahn- Teller lattice distortions that are present in this 
manganite cause a canting instability of its magnetic 
groundstate, driving a reconstruction of its charge or- 
dering from site-centered towards bond-centered. The 



resulting non-collinear magnetic ordering induces in 
Lai/2Ca!/ 2 Mn03 a purely electronic polarization of sev- 
eral ^(C/cm 2 , almost two orders of magnitude larger than 
the one of a typical multiferroic such as TbMn 2 Os [5]- 

We consider the half-doped manganite, 
Lai/2Ca 1 / 2 Mn03 (LCMO) in the experimentally 
observed antiferromagnetic (AFM) CE double-zigzag 
spin state JT2J, [13] (see FigJIJi), which is stable below 
Tjv ~ 155 K. Despite having a long history [13] IT4] . 
the experimental crystallographic and corresponding 
electronic structure of LCMO (as that of the closely 
related Pri^^Ca^MnOa, x ~ 0.5) is still debated. 
On one hand, a traditional checkerboard CO state 
has been proposed [12], given by the alternation of 
orbitally-ordered Mn 3+ (at the center of Jahn- Teller 
distorted octahedra) and Mn 4+ (in largely undistorted 
octahedra) with a mostly site-centered (SC) CO. On 
the other hand, a bond-centered (BC) model has been 
suggested [15] [16], based on the so-called Zener-polaron 
(ZP) state where equivalent Mn d A ions, showing no 
charge-disproportionation (CD), couple into ferromag- 
netic (FM) dimers sharing a spin-polarized hole on the 
intermediate O atom. In both SC and BC cases, we find 
magnetic groundstate structures that break inversion 
symmetry (IS) and result in polar states with relatively 
strong ferroelectricity [T7] , 

First-principles studies have proven to be tremen- 
dously helpful to shed light on the microscopic ori- 
gin and on the quantitative evaluation of the elec- 
tric polarization P in several improper magnetic ferro- 
electrics [HI HH |20j [21]. Our density functional the- 
ory (DFT) simulations are performed within the general- 
ized gradient approximation (GGA) [35] to the exchange- 
correlation potential and treating the Mn d electrons 
via a Hubbard-like potential within GGA+U [23] (un- 
less otherwise noted, U = 4 eV, J = 0.9 eV). We used 
the Vienna Ab-initio Simulation Package (VASP) [24] . 
including the non-collinear-spin formalism |25] and the 
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FIG. 1: (Color online) Spin-directions in the MnC>2 plane for 
a) 9 =0° (AFM-CE) and b) 9 = 90° (_L). Zig-zag chains and 
FM dimers highlighted, c) Electronic polarization P e u (in 
liC/Cm 2 ) vs9 (in de grees) for U = 4 eV and U = 8 eV for the 
experimental ionic structure of Ref. [T2| (CO) and of Ref. [27] 
(ZP). "Half-doping" is simulated both with i) a checkerboard 
arrangement of La and Ca ions and ii) by adding half extra- 
electron for each Ca ion in CaMnOa and compensating with 
a homogeneous positive background (denoted as N*) 



Berry-phase (BP) approach to evaluate electric polar- 
ization P [26]. The cut-off for the plane- wave basis set 
was chosen as 400 eV for the collinear and non-collinear 
spin-configurations and a [3,3,4] mesh was used for the 
Brillouin-zone sampling. In the BP approach, we inte- 
grated over 12 k-space strings parallel to either the a or 
b axis, each string divided in 8 k-points. The experimen- 
tal lattice parameters and ionic positions were taken from 
Ref. [12] for the CO-like structure (P2 1 /m space group) 
and from Ref. [57] for the ZP-like structure (P2\nm space 
group) . 

From our calculations we find that the CE-type AFM 
(Fig. [iji) is insulating (see Fig. [2^) and clearly or- 
bitally ordered, with a small CD, consistent with pre- 
vious collinear-spin electronic-structure works [25]. For 
U = 4 eV, there is a clear gap, the CD amounts to 8 ~ 
0.15 e~ whereas the magnetic moments are ~ 3.3 \xb 
and 3.05 (ig on the nominally Mn 3+ and Mn 4+ , respec- 
tively. The BP calculation of P shows, as expected, that 
the centrosymmetric CE-type structure is paraelectric. 
However, a spin-rotation immediately induces a ferro- 
electric moment. We consider a rotation of the spins 
of two neighboring Mn (one 3+ and one 4+) along the 
up spin-chain by an angle 9 and, correspondingly, one 
dimer in the down spin-chain by the same angle. This 
particular spin rotation is motivated by the fact that it 
tunes the magnetic CE state continuously towards the 
one compatible with ZP structure (denoted as _L in what 
follows), which corresponds to 9 = 90° and is shown in 
Fig. [T}d. The calculated electronic polarization P e i e in- 
creases monotonically with 9, reaching ~ 3 [iG'/cm 2 for 
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FIG. 2: (Color online) DFT results: a) Total DOS for differ- 
ent 9: DOS are arbitrarily shifted on the y-axis with a shift 
proportional to 9. The zero of the energy scale marks the 
Fermi level, b) Isosurface of the e g bands in the _L structure: 
view on the Mn02 plane. 

9 = 90°. Importantly, the results are stable with respect 
to an increase in the value of U, and to the details of how 
the half doping is achieved within DFT (see Fig. [lj). 

The resulting total density of states (DOS) shows that 
a larger value of 8 broadens the e g band ( 0.7 eV wide), 
therefore reducing the band-gap (see Fig. [2^). In Fig.[2]3 
we plot an isosurface of the e s -bands, indicating a clear 
orbital ordering (OO), with two kinds of different Mn: 
the nominal Mn 3+ shows a staggered (3x 2 —r 2 )/ (3y 2 — r 2 ) 
orbital arrangement, whereas the nominal Mn 4+ shows a 
much more isotropic charge-distribution, as given by par- 
tial occupation of both (3z 2 — r 2 ) and (x 2 — y 2 ) orbitals. 
Since this situation is by far similar to the collinear CE 
case (not shown), it proves that the spin-rotation alone 
does not alter significantly the CO/OO. 

In order to investigate the stability of the CE-type 
AFM state with respect to the rotation of spin-dimers, 
we study the degenerate double-exchange model in the 
presence of the inter-orbital Hubbard repulsion and Jahn- 
Teller (JT) lattice distortions, with the Hamiltonian, 

a/3 

H = -^2t^cos(e i3 /2)4 a c jl3 + Uj2n. la n lb 
<«) i 
+J s ^S l -S J -A^Q J .x. i + |^Q 2 , (1) 

(ij) 1 1 

where c (c^) is the annihilation (creation) operator, and 
a, (3 are summed over the two Mn-e g orbitals d x 2_ y 2(a) 
and d^ z 2_ r 2{b). t % £p denote the nearest-neighbor hop- 
ping amplitudes: t x aa = t y aa = t, tf b = t\ h = t/3, 
t x ab = tf a = -VV3, t\ = t\ a = t/V3. The 4J are the 
angle between neighboring Mn-t2 g spins Si and Sj and J s 
is their AFM superexchange. A denotes the strength of 
the JT coupling between the distortion = (Qi X ,Qiz) 
and the orbital pseudospin -rf = J2 a f3 c L-^q/3 c «/3' where 
r M are the Pauli matrices. K is a measure of the lattice 
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FIG. 3: (Color online) Model-hamiltonian results: a) Total 
energy E vs 8 for different values of A. The 9 dependence of 
b) DOS-gap A, and c) the bandwidth W of the occupied e 9 . 



stiffness (set to unity). Energies are in units of t (esti- 
mated ~ 0.2eV from the e g bandwidth in Fig. [2| 



In 

the absence of U, a CE state with CO/OO is found to 
be the groundstate over a wide regime in the parameter 
space [S3]. In the absence of A (and for U > 10), we find 
that the _L state has lower energy than the CE state (see 
Fig. [3^) [30]. Inclusion of the JT coupling leads to the 
stability of a state with an intermediate 9. Eventually, 
beyond a ([/-dependent) critical value of A, one recovers 
the CE state as the groundstate of this model, due to the 
more pronounced energy gain via JT distortions in the 
CE state as compared to that in the _L state. 

It is interesting to note that when computing the 
spin-structure factor for various 9, we find that there 
is no change in the location or strength of the major 
peaks, positioned at (it, 0) , (0, tt) and (tt/2, 7r/2). This 
makes it very hard to differentiate between the magnetic 
structures corresponding to different 9 from simple neu- 
tron scattering experiments. We also find a decrease 
in the DOS gap (Fig. ^jp), and a broadening of the e g 
bands (Fig. |3p) upon increasing 9, similar to the first- 
principles results (see Fig. [2^). This can be understood 
in terms of the effective hoppings arising via the DE fac- 
tor cos(Ojj/2). In the CE state each chain has a perfect 
hopping, but there is no inter-chain hopping. However, 
the spin-dimerized state allows for inter-chain hopping 
at the cost of reduced hopping within chains, leading to 
an overall gain by a factor (1 + v / 2)/2, which very well 
explains the increase in the bandwidth W (Fig. [3J;) . 

Let us now focus on the multiferroicity. The creation 
of non-collincar "spin-dimcrs" along with the small but 
detectable CD breaks IS (present in the CE spin ar- 
rangement), allowing FE polarization along the a-axis 
as a realization of the intermediate BC/SC CO-picture 
proposed by Efremov et al. [10 . There is however one 
important difference between our work and those model 
predictions |10j : there, the spin-rotation is sufficient to 
progressively transform a SC-CO into a BC-CO (ZP) 
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FIG. 4: (Color online) a) DFT isosurface of the e g bands in 
the LT-0 structure. ZP are highlighted, b) Atomic displace- 
ments (not in scale) of the LT-0 Mn02 plane. Atoms marked 
with the same color are structurally equivalent. Zig-zag spin 
chains highlighted. 



and, as such, the _L-structure regains centrosymmetry 
(i.e. leading to the expectation of P e i e = 0). In our case, 
the spin-rotation - even for the largest 9 = 90° - is not 
a strong enough factor to produce a BC situation, due 
to the structural inequivalency between the two kinds 
of Mn which governs their charge distribution. Indeed, 
in our DFT simulations, the larger the deviation from 
the collinear CE-type, the stronger is the "asymmetry" 
introduced by the spin-dimers along the chain and the 
closer one gets to the ideal realization of the interme- 
diate BC/SC situation, explaining why the maximum is 
located at 9 = 90°. According to this mechanism, CD 
is an essential ingredient in the rising of P: there would 
be no ferroelectricity without CD (despite Mn cations 
keeping a centrosymmetric distribution). Rather, it is 
the combination of spin-rotation and CO that induces 
ferroelectricity. 

We now discuss an alternative lattice structure in 
which the atomic arrangement is characterized by a 
structural Mn-Mn dimerization. In particular, Rodriguez 
et al. [27 j proposed two different LCMO structures (re- 
ferred in that paper as LT-M and LT-O): the LT-M shows 
basically the same configuration and symmetries as in 
Ref. [IS], whereas the LT-0 shows neighboring octahe- 
dra in which both Mn are off-centered and with "long" 
MnO bonds directed along the same Mn-O-Mn line. This 
is characteristic of a ZP-like structure [15], in which the 
two Mn (despite being still inequivalent by symmetry) 
are expected to appear electronically more similar than 
in the SC CO/OO model. Indeed, when plotting the 
LT-0 isosurface for the e g manifold (see Fig. |4^i), the 
two kinds of Mn show a very similar charge distribu- 
tion, as also confirmed by the small difference between 
their moments (about 3.2 /j,b and 3.32 (J-b)- The peculiar 
OO shows pairs of Mn with 3x 2 — r 2 orbitals alternated 
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with pairs of Mn with 3y 2 — r 2 orbitals, in agreement 
with the ZP picture. However, from our calculations 
we do not find an appreciable spin-polarization on the 
O in the ZP-like dimer, which is ~ 0.05/is, at variance 
with much larger values predicted by Hartree-Fock cal- 
culations [3TJ [35]. Therefore, small charge-transfer ef- 
fects are expected, calling for further studies focused on 
the O spin-polarization (especially in Pri-^Ca^MnOs, 
0.3 < x < 0.5, where the ZP seems the ground state). 
In this regard, we note that with AFM-CE spin config- 
uration the LT-M shows a lower total energy (by ~ 30 
meV/Mn) than the LT-O, suggesting as the ground state 
for LCMO a mostly SC state with CD. 

The calculation of the FE polarization for this struc- 
ture in the AFM-CE spin configuration gives particu- 
larly interesting results. In the P1\nm space-group, 
a non-switchable polarization is allowed along the a- 
axis [IB]. However, we here focus on the 6-axis, investi- 
gating the possibility of magnetically switchablc multifcr- 
roicity. The 2i symmetry, element of the P1\nm space- 
group, forbids any ionic component of P along the 6-axis: 
indeed we find P\ on =0. However, when calculating the 
electronic contribution, a large value is obtained (P% e = 
7.2 /xC/cm 2 ), which should be easily experimentally de- 
tected in untwinned high-quality crystals. To investigate 
the origin of P, we show in Fig. |4Jd the displacements of 
the atoms with respect to a reference centrosymmetric 
Puma structure [12] . It is clear that the components of 
the displacements along the a-axis do not cancel, result- 
ing in a finite P a . However, along b, there is a complete 
cancellation in the displacements, consistent with the 2i 
symmetry, with the ionic P summing up to zero. Never- 
theless, when imposing the AFM-CE spin configuration, 
the 2i screw axis is not any more a symmetry opera- 
tion in the magnetic space-group. For example, the O 
atoms labeled as On,i, are alternatively bonded to Mn 
with parallel (Ofj 4 ) and with antiparallel (Ojj 4 ) spins: 
they are therefore structurally equivalent but electroni- 
cally incquivalent, as suggested by the e g charge density 
plot (Fig. [iji). To further substantiate the inequivalency 
of the O as the source for Pj! le , we created spin-dimers 
by coherent rotations. Indeed we find that P decreases 
continuosly upon increasing 9 and vanishing at 90° (see 
Fig.nl as the inequivalence of O v n 4 and Ojj 4 disappears. 

In summary, we have found two different mecha- 
nisms for improper ferroelectricity in LCMO. In the first 
one, starting from centrosymmetric ionic positions and 
double-zig-zag ferromagnetic spin-chains, an intermedi- 
ate bond-centered/site-centered polar charge distribution 
is achieved by means of a spin-dimerization. This in- 
duces a dramatic ferroelectric response with P up to few 
^tC/cm 2 . The second mechanism, active in the bond- 
centered ZP-like lattice structure, induces ferroelectricity 
along 6-axis, with the AFM-CE spin arrangement lifting 
the 2i symmetry and paving way to a value of Pb which 
is largest to date in the whole class of improper magnetic 



ferroelectrics. Recent electric field gradient experiments 
on doped manganites have also suggested the existence of 
ferroelectric domains in the charge ordered regime |33j . 
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